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Abstract. Doubly non-negative matrices arise naturally in many setting including Markov random 
fields (positively banded graphical models) and in the convergence analysis of Markov chains. In 
this short note, we settle a recent conjecture by C.R. Johnson et al. [Linear Algebra Appl. 435 
(2011)] by proving that the critical exponent beyond which all continuous conventional powers of 
n-by-n doubly nonnegative matrices are doubly nonnegative is exactly n — 2. We show that the 
conjecture follows immediately by applying a general characterization from the literature. We prove 
a stronger form of the conjecture by classifying all powers preserving doubly nonnegative matrices, 
and proceed to generalize the conjecture for broad classes of functions. We also provide different 
approaches for settling the original conjecture. 



1. Introduction and main result 

The study of operations preserving different notions of positivity is an important topic in matrix 
analysis (see e.g. Bhatia [Ij, Horn and Johnson |3]). The purpose of this short note is to prove the 
critical exponent conjecture for doubly nonnegative matrices formulated recently by Johnson et al. 
[5]. A real symmetric matrix A is said to be doubly nonnegative if it is both positive semidefinite, 
and entrywise nonnegative. Given a real function / : [0, oo) — >• [0, oo) and a doubly nonnegative 
matrix A, we define the matrix f{A) using the spectral decomposition of A. Namely, if ^ = UDU^, 
where U is orthogonal and D = diag((iii, . . . , (i„,„) is diagonal, define f{A) := Uf{D)U'^, where 
f{D) := diag(/(dii), . . . ,f{dnn))- Assuming A is doubly nonnegative, it is very natural to seek 
conditions under which f{A) is also doubly nonnegative. 

A similar question can be asked when the function / is applied entrywise to the elements of A. 
Denote by f[A] := {f{aij)). In the particular case where f{x) = x" for < a G M, it has been 
shown by FitzGerald and Horn [3| (see also Horn and Johnson [H Chapter 6.3]) that f[A] is doubly 
nonnegative for every doubly nonnegative matrix A if and only if a GNU[n — 2,oo). Note that in 
each of the two questions above, one of the two nonnegativity conditions is merely used to ensure 
that f{A) or f[A] is well-defined, and the real question is whether the other notion of nonnegativity 
is preserved. 

Johnson et al. are interested in the first question above for the power functions f{x) = x" for 
real a. Namely, which conventional powers preserve the set of doubly nonnegative matrices of a 
given order? The following phase transition can be shown to occur for every integer n > 0. 

Theorem 1.1 ([21 Theorem 2.1]). There is a function /i : N ^ [0, oo) such that for any n-by-n 
doubly nonnegative matrix A, the matrix f{A) = A" is doubly nonnegative for a > n{n). 
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The smallest such /i(n) always exists by continuity and is denoted by m{n). The quantity m(n) is 
called the (conventional) critical exponent (see [5j), and it is of great interest to identify it. Johnson 
et al. [5j provide a lower bound for m(n) for all n, and a sharp upper bound for small values of n. 

Theorem 1.2 ([5], Theorems 3.1, 3.2, 4.1]). The critical exponent m{n) > n — 2. If n < 6, then 
the critical exponent m{n) = n — 2. 

The authors then proceed to conjecture if the above result is true in general: 

Conjecture 1.3 (O Johnson, Lins, and Walch, 2011]). For all integers n > 2, the critical exponent 
m{n) = n — 2. 

The goal of this note is to settle this conjecture and prove that indeed m(n) = n — 2 for all 
1 < n E N. More generally, we prove this conjecture (in fact, a stronger statement of it) for a larger 
class of functions. In order to state our main result, we first need some notation. 

Definition 1.4. Let n > 2 be an integer. Denote the set of doubly nonnegative n x n matrices by 
DNn- Now given a function / : [0, oo) [0, oo), define the critical exponent for / to be 

(1.5) m{f,n) := inf {0 < t e M : ^"/(^) £ DNn € DNn,ya > t}. 

We now state the main result of this paper. 

Theorem 1.6 (Critical exponent conjecture, stronger form). Given u > and /3 € M, the critical 
exponent for f (x) := {x + u)~^ ism{f,n) = n — 2 + max(/3, 0). Moreover, below the critical exponent 
there are only finitely many values < a < m{f,n) such that the map g[x) = x°^{x + u)~^ preserves 
DNn, diT'd these are contained in the set 

(1.7) {m + max(/3, 0) : m = 0, 1, . . . , n — 3}. 

Remark 1.8. Note that the critical exponent conjecture 11.31 for powers by Johnson et al. [5] is 
the special case where /3 = 0. Additionally, Theorem 11.61 strengthens the statement of the critical 
exponent conjecture significantly in two ways: (1) it discusses the behavior of A^'f^A) for values of 
a lower than the critical exponent, and (2) it considers more general functions than conventional 
matrix powers. 

We will also see in Section [2.11 that showing that g{x) = x^~^{x + u)~^ preserves DNn for all 
u > actually implies the critical exponent conjecture. This reason also motivates our trying to 
prove Theorem 11.61 which is a more general result than the standard form of the critical exponent 
conjecture. 

2. Proof of the Strong critical exponent conjecture 

We now settle Conjecture 11.31 bv proving the more general result given by Theorem 11.61 To do 
so, we shall employ the following result from the literature. 

Theorem 2.1 ([6l Micchelli and Willoughby, Corollary 3.1]). Suppose g{x) > is continuous for 
X > 0. Then g preserves DNn if o-'^d only if all the divided differences of g in [0, oo) of order 
1, 2, . . . , n — 1 are nonnegative. 

Theorem 12.11 is shown by replacing g by the Newton interpolation polynomial taking values g{\i) 
at the eigenvalues < Ai < • • • < A„ of A, and proving that each summand is doubly nonnegative, 
i.e., 

(2.2) {A - All) ■■■{A- Afcl) G DNn ^A G DNn, l<k<n-l. 

Note that if g is C"""*^, then g preserves DNn if and only if g^'^\x) > for all < i < n and 
X > (by the mean value theorem for divided differences). As a consequence, we have the following 
corollary. 
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Corollary 2.3. Conjecture \1.3\ holds for all integers n > 2. Moreover, g{x) = preserves DNn 
if and only if a £ N U [n — 2, oo). 

Proof. Immediate by considering g{x) = x". □ 

Note that the set of powers which preserves DNn for conventional powers (Corollary I2.3P is exactly 
the same as for Hadamard powers p]: a G N U [n — 2, oo). 

Remark 2.4. Besides settling the critical exponent conjecture for powers, Theorem 12.11 can also 
be used to establish the existence of a critical exponent in very general settings. For instance, a 
family of smooth functions {fa)a>o will have a "critical index" m{n) for DNn, if and only if 

(2.5) fi'Kx)>0, Vx>0, « = 0,l,...,n-l 

for all a large enough. A stronger statement would be that (j2.5p holds if and only if a G S{n) U 
[m(n), oo) for some finite set S{n). We show in Theorems 11.61 and 13.11 that this is indeed the case 
for large families of functions. 

To prove Theorem 11.61 using Theorem 12. H we need the following preliminary result. 

Lemma 2.6. For all a, /3 € M, u, x > 0, and < r G Z, 

dx'-ix + u)/^ (x + u)P+'-^\iJ -^^11^ ^ 



r\ r(a + l)r(a - f3 - i + 1) 



(x + u)P+'- ^\i T(a-i + l)T(a - /3 - r + 1) ' 

Proof. The proof is by induction on r > 0; the base case is clear. Now given the result for r, 
differentiate the right-hand side using the quotient rule. A typical summand involves differentiating 
{x+'u)0+- ' derivative is 

(x + ti)^+''(a - i)x"-*-i - (/? + r)(x + u)^+'^-ix°-* 



^a—i—l 



(X + u)2{/3+0 

^ ((a - i){x + u) - {/3 + r)x) 



{x + u 

^ {{a - (3 - r - i)x''-' + (a - i)ux'^-'-^' 



(X + 7/)/3+^'+l 

Now multiply this quantity by the coefficient (J'-^u'^Tr^i, where 

r(a + l)r(a-/3-i + l 



' T{a-i + l)T{a- (3-r + l)' 

Add the products from z = to r, and collect like powers of u to get the (r + l)th derivative. 
We claim that this is indeed of the claimed form, which would complete the induction step and 
show the result. But this is a straightforward verification for each term: the coefficient of in the 
(r + l)th derivative comes from only the i = term: 

x"(a — /S — r) _ x" 
as desired. Similarly, the coefficient of u'^'^^ also comes from only the i = r term: 



4 



DOMINIQUE GUILLOT, APOORVA KHARE, AND BALA RAJARATNAM 



It remains to show that the coefficient of -3- — rr (for 1 < z < r) in the sum in the (r + l)th 

derivative, which equals {^~^^^Tr-\-i^i, comes from adding up the two terms in the above sum of 
derivatives. In other words, it suffices to show that 



) (a - /3 - r - i)Tr,i + ( / ^) (« " ^ + ^Wr,i^i = ^ ^) ^r+l,i. 



Now note that each of the three terms (^.^-^Tr^i-i, (^Fr^i, [^~^^^Tr-\-i^i have the following expression 
common: 

^ ' j=0 j=i 

Thus, it suffices to show the equality after having "cancelled" this common expression from all 
terms - namely, to show that 

{r - i + l){a - P - r - i){a - i + 1) + i{a - i + l){a - 13 - i + 1) 
= (r + l)(a - i + l)(a - /3 - r). 

Once again, it suffices to show this without the common expression (a — i + l) in each of the terms. 
But this is a straightforward computation. □ 

Finally, we now prove the main result of this note. 

Proof of Theorem \1.6\ Suppose first that a > n — 2 + max(/3, 0). Then for all < i < n and x > 0, 
the ith derivative 



dx^ {x + u)p 

is indeed defined and can be computed using Lemma 12.61 It is easy to check that all of these 
derivatives are nonnegative, since x,u > 0, and a — {n — 2) and a — /? — (n — 2) are nonnegative. By 
Theorem 12.11 (see the remarks preceding Corollarv 12. 3p . we infer that x'^{x + u)~^ preserves DNn 
for all a > n — 2 + max(/3, 0). 

To complete the proof, suppose now that 0<a<n — 2 + max(/3, 0) and a is not in the set (jl.7p . 
Define g{x) = x'^{x + u)~^ : [0, 00) — )• [0, 00). By Theorem 12. H it suffices to show that there is at 
least one < i < n and x > such that g^^\x) < 0. There are two cases: 

(1) Suppose (3 < 0. Now compute by Lemma [261 



a—m—2 



{x + n)/5+™+2 



m+l , 1 9\ m+1 

j=0 i<m+2 ^ ^ j=0 j=i 



Now if Q € (m, m + l) for some integer m € [0, n — 3], then u™"'"^ Y\^=Qi^~3) smallest 
degree term in x and is negative. Therefore g^^'^'^\x) < for small x > 0. 
(2) Suppose /3 > 0. Once again using Lemma [2^ 



a—m—2 



(x + uy+' 



m+2 



m+l / I 9\ m+l 

X-+2 n(«-/3-j)+xj;r+ v^^^'~' n (" - n ( 

i=0 i>0 ^ ^ ^ j=0 j=i 



Now if a — /3 € (m, m+l) for some integer m G [0, n — 3], then a similar analysis as in the 
previous case, but this time with x 00, implies that g^^~^'^\x) < for some x > 0. 

□ 
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2.1. Alternate proof. We now provide an alternate proof of the original critical exponent conjec- 
ture in [5j. Similar to Johnson et al.'s result in [5] (see Theorem II ■2p . we only show the conjecture 
(i.e., that m(/, n) = n — 2 for f{x) = 1) for small values of n, although we also indicate an approach 
for all n. 

The first simplification is in showing that to verify whether or not a given function / preserves 
DNn, it is enough to focus on a particular off-diagonal entry: 

Lemma 2.7. Given f : [0, oo) — > [0, oo) and a fixed pair of integers 1 < i ^ j < n, the matrix f{A) 
is doubly nonnegative for all A € DNn if o-nd only {f{A))ij > for each A € DNn- 

Proof. One implication is immediate from the definition. Conversely, given any A = UDU^ with 
U unitary and D diagonal, note that f{A) = Uf{D)U^ is positive semidefinite. Hence its diagonal 
entries are nonnegative. We now consider the off-diagonal entries of f{A). Note that for every 
n-by-n permutation matrix P, 

fiPAP^) = f{{PU)DiPUf) = {PU)f{D){PUf = Pf{A)P^. 

Since conjugating by P preserves DNn, hence {Pf{A)P^)ij > for all P,A. This is equivalent to 
saying that f{A)i/j/ > for all i' ^ j' ■ □ 

Next, note by Theorem 12.11 and Lemma 12.61 that for every n G N, 
(2.8) A'^-^iA + uI)-^ e DNn, yA£ DNn,u>0. 

We now show that this fact alone implies the critical exponent conjecture 11.31 

Proposition 2.9. Fix an integer n = hq >2. Then Equation (j2.8|) for n = uq implies the critical 
exponent conjecture \1.3\ for n = uq. 

Proof. By [5, Theorem 3.2], we know that m(no) > — 2. To show that m{no) = uq — 2, assume 
first that A € DNn^ is nonsingular. Clearly, A^ is doubly nonnegative if A € DNn^ and A; € N. 
Therefore, it suffices to prove that A'^ is doubly nonnegative for every k < q < k + \ and every 
integer k > uq — 2. To show this, note that for k < q < k + 1 and x > the following formula 
holds: 

(2.10) x^ = '^"^^^"^^""^ r f-'^-^x'^+^x + t)-' dt. 

71" Jo 

Equation (j2.10p is shown by a standard contour integration followed by an application of the 
residues theorem (see e.g. [21 Chapter V, Example 2.12]). But now the matrix = UD'^U^ can 
be written as: 

(2.11) A'^ = r t^-^-^A'^\A + ti)-' dt = ^M^^m r ,(^) dt, 

Jo Jo ' 

where fk,t{x) ■= x'^^^{x + t)~^ for t > and x > 0. Note that sin{{q — k)TT) > since < q — k < 1. 
Moreover, since for every A € DNng and k > hq — 2, the matrix fk^t^A) is doubly nonnegative by 
Equation (|2.8p . it follows immediately that A'^ is doubly nonnegative. 

Finally, if A is singular, then A*^ = lim^_^Q+{A + ul)'^ and the result follows. □ 

Remark 2.12. Equation (j2.1ip is frequently used in the study of positive definite matrices. See 
e.g. Bhatia [H Theorem 1.5.8]. 

Given Proposition 12.91 it is thus of interest to prove Equation (j2.8p from first principles in order 
to provide an elementary proof of the critical exponent conjecture in [5j. We conclude this section 
with a case-by-case verification of Equation (|2.8p for n < 3. Such an approach (but using different 
arguments) was adopted by Johnson et al. in [5] when they showed the conjecture for n < 6. 

Proposition 2.13. The critical exponent conjecture \1.3\ holds for n < 3. 



6 



DOMINIQUE GUILLOT, APOORVA KHARE, AND BALA RAJARATNAM 



Proof. In light of Proposition 12.91 it suffices to show that Equation ()2.8p holds (for n < 3). First 
note that A^~^{A + ul)~^ is clearly positive semidefinite. By using Lemma 12.71 it now suffices to 
show that the (l,2)-entry of A'^~^ {A + ul)~^ is nonnegative. Denote this entry by ai2,n- It is then 
easily shown that for A G DN2, 

ai2U 

= det{A + ul2) - 

Similarly, for n = 3, 

ai2,3 = det{A + {a^ det A + ai2u{Cu + C22 + ^33) + u^{ai2an + 012022 + 013023)) , 

where Cu are the principal cofactors/minors of A, hence are nonnegative. This concludes the proof 
since all entries of A are nonnegative. □ 

3. Generalization of the strong critical exponent conjecture 

Note that Theorem 12.11 can be used to systematically study and resolve the (strong) critical 
exponent conjecture for broad classes of functions. We conclude this note by showing such a result. 

Theorem 3.1. Let f : [0, 00) [0, 00) be a C"'^^ function satisfying 

(3.2) /(x)>0, lim /(x) > 0, f^'\x)>0, lim /W(x) = Af^ < 00, Vx > 0, < i < n. 

>0+ X— >0+ 

Then the values < a € M such that x"f{x) preserves DNn are N U [n — 2, 00). In particular, the 
critical exponent for all such functions f is m{f, n) = n — 2. 

Proof. Set g{x) := x"f{x). In order to use Theorem 12. H we evaluate g^'^\x) using the product rule: 

i-i 



^7«(x)=EC)n("-^>""'/^ 



^Xj 

If a S N then every product nj=o('^ ~ j) that includes a negative term also includes 0, whence all 
derivatives of g are nonnegative at x > 0. Similarly, if a > n — 2, then g^^\x) > for x > 0. Finally 
if a ^ N and a < n — 2, assume a € (m, m + 1) for some integer m E [0, n — 3]. Then the I = m + 2 
term for g^'^^'^\x) is negative, and hence g^"^~^'^\x) < for small values of x > by ()3.2p . □ 

Remark 3.3. Note that the critical exponent conjecture in [5J is again a special case of Theorem 
13.11 - with /(x) = 1. Moreover, the cone ^ of functions / satisfying (13. 2p is a very large family that 
contains all absolutely monotonic functions with nonzero constant term as well as exponential-type 
functions such as e^' and cosh(x). Additionally, ^ is closed under multiplication, exponentiation, 
composition, and taking nontrivial nonnegative linear combinations. 
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